Coherent intense resonant laser pulses lead to interference in the time domain 
observable in the spectrum of the emitted particles 
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The dynamics of atomic levels resonantly coupled by a coherent and intense short high-frequency 
laser pulse is discussed and it is advocated that this dynamics is sensitively probed by measuring the 
spectra of the particles emitted. It is demonstrated that the time-envelope of this laser pulse gives 
rise to two waves emitted with a time delay with respect to each other at the rising and falling sides 
of the pulse, which interfere in the time domain. By computing numerically and analyzing explicitly 
analytically a show-case example of sequential two-photon ionization of an atom by resonant laser 
pulses, we argue that this dynamic interference should be a general phenomenon in the spectroscopy 
of strong laser fields. The emitted particles do not have to be photoelectrons. Our results allow 
also to interpret the already studied resonant Auger effect of an atom by intense free electron laser 
pulses, and also to envisage experiments in which photons are emitted. 

PACS numbers: 33.20.Xx, 41.60.Cr, 82.50. Kx 



The interaction of an atom with intense laser fields 
has been widely studied. If the field is essentially 
monochromatic, the physics is well described by a time- 
independent Hamiltonian in the basis of 'dressed' elec- 
tronic states or Floquet states (see, e.g., Refs. (H-Q). The 
inclusion of relaxation mechanisms, such as autoioniza- 
tion or subsequent ionization, gives a 'dressed' state a 
finite width, and it becomes unstable 5]. The concept 
of dressed states is applied in practically every branch 
of spectroscopy of optical lasers operating in the nano 
and picosecond regimes. If the laser pulses are shorter, a 
Floquet basis is still useful, but one has to take the time- 
dependence of the pulse explicitly into account. Many 
new phenomena arise due to the impact of this time- 
dependence 

One class of such phenomena extends the well-known 
stationary Rabi-doublets existing in strong fields owing 
to ac-Stark splitting or Autler-Townes effect [11]. Be- 
cause of the short optical pulse, the value of this split- 
ting varies, resulting in the appearance of a multiple- 
peak interference pattern in the computed autoionization 
[12| and resonant multiphoton ionization 1^, 14 1 electron 
spectra. This new pattern is attributed to the temporal 
coherence of a pulse strong enough to induce Rabi oscil- 
lations between resonantly coupled states 12Ul4j. How- 
ever, a physically simple explanation of the phenomenon 
is still missing [15] . 

Although these theoretical predictions are of relevance 
and were made a long time ago, they have not been 
verified experimentally so far. To our opinion, this is 
due to the optical regime. First, in this regime there 
are rarely well separated resonances and there is often a 
dense spectrum of close-by Rydberg and doubly-excited 
states which also participate in the dynamics. Second, 
these states induce additional ac-Stark shifts which vary 



in time [lfj . Third, one is often in the vicinity of ion- 
ization thresholds and ionization is particularly efficient 
there. All of these additional states and effects strongly 
smear out the pronounced interference pattern which 
would be obtained if only two or three states were reso- 
nantly coupled by the pulse. 

The situation becomes particularly promising by the 
advent of the new generation of light sources, like attosec- 
ond lasers 0], high-order harmoni c g eneration sources 
17 , lij , and free electron lasers 19|, 20] to produce ultra- 
short and intense coherent laser pulses of high frequen- 
cies. The above mentioned shortcomings which impede 
experimental verifications by optical pulses, are absent 
at higher frequencies and one can study the dynamics of 
a few well separated electronic states (e.g., core-excited 
states) resonantly coupled by a short coherent pulse. Un- 
less the intensity is very high, the resonant dynamics will 
not be affected by ac-Stark shifts arising from nonessen- 
tial states, and the impact of direct ionizations is not 
substantial since the photoionization probability usually 
decreases with the photon energy. We thus concentrate 
in this work on the high-frequency regime and discuss a 
fundamental consequence of the nature of intense coher- 
ent laser pulses on spectroscopic observables. Due to the 
high carrier frequencies, much of the physics follows the 
evolution provided by the pulse envelope nearly adiabat- 
ically up to rather short pulse durations. This makes the 
underlying physics particularly transparent. 

Let us consider two bound electronic states of an atom 
coupled resonantly by a strong laser pulse (the carrier 
frequency can be different from the field-free resonant 
frequency to compensate for the emerging energy detun- 
ing by the AC Stark effect). The two initially degenerate 
'dressed' states repel each other by the field-induced cou- 
pling and split in energy. If the pulse envelope supports 
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many optical cycles of high frequency, the field-induced 
coupling between the two electronic states adiabatically 
follows the pulse envelope [l6| . Consequently, the energy 
splitting will adiabatically increase when the pulse ar- 
rives and then decrease when the pulse expires. If the 
atom emits particles during its exposure to the pulse 
(photoelectrons, Auger electrons, photons), it will be- 
come evident below that the particles emitted when the 
pulse rises have the same kinetic energy as those emit- 
ted when the pulse decreases. The respective two waves 
emitted with a time delay with respect to each other will 
interfere and their spectrum will exhibit a pronounced 
interference pattern. We would like to call this kind of 
interference, dynamic interference. 

Although we concentrate here on high-frequency short 
pulses coupling two bound states, we mention that 
bound-continuum coupling by such pulses also leads to 
dynamic interference in the ionization spectra of atoms 



2l| and model anions 22j- Furthermore, oscillations in 



the total multiphoton ionization yield as function of laser 
intensity have been observed for atoms exposed to opti- 
cal lasers 23], l24[ and interpreted as arising from interfer- 
ences of electrons emitted at different times 24j . We shall 
demonstrate here that dynamic interference is a general 
consequence of the finite nature of intense high-frequency 
laser pulses, and leads to pronounced patterns observable 
in the spectrum of the emitted particles. We first con- 
centrate on a show-case example of sequential two-photon 
ionization of an atom by strong pulses. The example is 
of much interest by itself, since the coupled two-level sys- 
tem is probed here by a second photon of the same pump 
pulse. Our results pave the way for experiments on dy- 
namic interference by available laser pulse sources. We 
also briefly discuss the effect of dynamic interference in 
other branches of laser spectroscopy. 

Below, we consider an atom initially in its ground elec- 
tronic state 1 7) of energy of E] = chosen as the origin 
of the energy scale, which is resonantly excited into the 
intermediate state \R) of energy E R by absorption of a 
single photon and subsequently ionized by a second pho- 
ton into a final electron continuum state \Fe) of energy 
IP + e. Here, IP = Ep — Ej is the ionization poten- 
tial and s is the kinetic energy of the emitted photoelec- 
tron. Employing a linearly polarized coherent laser pulse 
£(t) = £q g(t) cosuit with pulse envelope g(t), the total 



wave function function of time reads [2lL I25W2 



= ai(t)\I) + a R (t)e- iut \R) + / a £ (t)e~ 2lu ' t \Fe)de 



where aj(i), a R (t), and a E (t) are the time-dependent am- 
plitudes for the population of the \I), \R), and \Fe) lev- 
els, respectively. The stationary states \R) and \Fs) have 
been 'dressed' by multiplying with the phase factors e ttJ}t 
and e 2lut [26] , to simplify the equations of motion. 
Inserting ^(t) into the time-dependent Schrddinger 
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FIG. 1: (Color online) Sequential two-photon ionization of H 
by a Gaussian-shaped pulse of 30 fs duration and resonant 
carrier frequency of u) = 3/8 a.u. = 10. 20 eV, which fits to the 
energy of the H(ls)-H(2p) excitation. Shown are the popu- 
lations of the ground state H(ls) and of the resonant state 
H(2p) as functions of the peak intensity after the laser pulse 
has expired. The vertical lines indicate the peak intensities 
at which the spectra depicted in Fig. [2] are computed. 



equation for the total Hamiltonian. and implying also 
the rotating wave 0, and local 2(| 3(1 31 1 approxi- 
mations, we obtain the following set of equations for the 
amplitudes (atomic units are used throughout) 



£>tg 

2 

-{E R - 



g(t)a R (t), 

^Tg 2 (t) -u)a R (t), 



idj(t) 

ia R (t) = -y5(t)oxW- 
io. e {t) = -g-atfC*) 

Here, D — (R\z\I) and d = (Fe\z\R) are the dipole tran- 
sition matrix elements for the excitation of the intermedi- 
ate state and for its subsequent ionization, respectively. 



(IP + £ - 2w) a e (t). 



(2a) 
(2b) 
(2c) 



The term — ^Tg 2 (t) in Eq. (|2bl) is the time-dependent 
ionization rate of the intermediate state responsible for 
the leakage of its population by the ionization into all 
final continuum states \Fe), turnin g th is state into a res- 
onance. Explicitly, V = 2Tr\d£ /2\ 2 [26|, |32|. 

To exemplify the present theory, we study the sequen- 
tial two-photon ionization of the hydrogen atom. In the 
process, H(ls) is resonantly excited to H(2p) state, which 
is then ionized. The photon energy was set to fit the exci- 
tation energy uj = E R = 3/8 a.u. = 10.20 eV. The com- 
puted dipole transition matrix elements for the excitation 
and ionization are D — 0.744 a.u. and d = 0.377 a.u., re- 
spectively. The system of Eqs. © was solved numerically 
employing a Gaussian pulse g(t) = e _t / T of r = 30 fs 
duration. Fig. Q] shows the populations of the ground 
state H(ls) and of the resonant state H(2p) after the 
laser pulse has expired as function of the peak intensity 
Iq = £q/87to!. The populations exhibit pronounced Rabi 
oscillations. To be noticed is that at the highest intensity 
considered in Fig. [T] the total photoelectron yield reaches 
just 7% indicating that the ionization by the second pho- 
ton is far from saturation. 

We now turn to the photoelectron spectra. For the cal- 
culations we have chosen the peak intensities at the max- 
ima of the ground state population indicated in Fig. [T]by 
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FIG. 2: Sequential two-photon ionization of H by a Gaussian- 
shaped pulse of 30 fs duration and resonant carrier frequency 
of uj = Er = 3/8 a.u. = 10.20 eV, which fits to the energy 
of the H(ls)-H(2p) excitation. Shown are the photoelectron 
spectra computed via Eqs. Q for different peak intensities 
indicated in the figure near each curve. The central electron 
energy so = 2ui — IP — 0.25 a.u. = 6.80 eV at which the 
photoelectron spectrum has its maximum in the weak-field 
case is indicated by a vertical line. 



vertical lines. At these intensities the atom manages to 
complete an integer number of Rabi cycles during the 
pulse duration. The spectra computed via Eqs. @ arc 
shown in Fig. [3] The spectrum computed for the low- 
est considered intensity of 5.2 x 10 11 W/cm 2 is rather 
close to that expected in the weak-field case, i.e., a Gaus- 
sian curve centered around £n = 2a; — IP = 0.25 a.u. = 
6.80 eV. As the field intensity increases and the atom 
manages to complete two Rabi cycles while the pulse 
is on (second spectrum from the bottom), the spectral 
distribution bifurcates, and is now minimal at £n. At 
the intensity 4.7 x 10 12 W/cm 2 when the atom has com- 
pleted three Rabi cycles, the spectrum bifurcates again 
and possesses now three maxima (third spectrum from 
the bottom). As the pulse intensity grows further, the 
spectrum continues to bifurcate again and again, exhibit- 
ing thereby distinct multiple-peak structures. Below we 
identify dynamic interference as the physical origin of 
these patterns. 

To start the discussion, we notice that the resonantly 
(lu = Eji~Ej) coupled dynamics of the |/) and \R) states 
in Eqs. (|2a[) and (|2b[) is governed by the 2x2 Hamiltonian 



H(t) 



At g(t) 
Ag(t) -|r ff 2 (i) 



(3) 



where A = ^f- . The ionization of the intermediate state 
by a second photon from the same pulse is described in 
Eq. © by the — ^Tg 2 (t) term, and actually probes this 



Hamiltonian. We may now follow the time evolution of 
the eigenvalues and eigenvectors of this Hamiltonian. 

When the pulse is on, the solution of Eq. © yields 
two decoupled resonances, which are superpositions of 
the initial \I) and intermediate \R) states: 



E ± (t)c±Ag(t)--Tg 2 (t), |±) 



ID. 

V2 



± 



V2" 



(4) 



This result is well justified if the pulse is not too strong 
and the ionization is far from saturation, i.e., when 
Ag(t) ^> \Yg 2 {t). These solutions describe two decou- 
pled time-independent resonances with time-dependent 
energies ±A g(t) induced by the field. Importantly, their 
energies move apart as the pulse arrives, and then move 
towards each other as the pulse expires. Both resonances 
are subject to the same leakage — jTg 2 (t), populating 
thereby the continuum states \Fe) via the ionization by 
a second photon. 

The decoupled resonances scenario enables one to un- 
cover the origin of oscillations in the spectra in Fig. [2] 
Using Eqs. Q we can rewrite the original Eqs. ((2]) in 
terms of the decoupled resonances |+) and |— ) and ob- 
tain the equations for the amplitudes a + (t) and a_ (i) of 
these resonances which can be solved analytically. Em- 
ploying the initial conditions a±(— oo) = 1/V2 we find 



a±(t) 



1 [TiAF(t)- 

y/2 C 



r/4 j(t)] 



(5) 



where F(t) = J^g^dt' and J(t) = f_ oo g 2 {t')dt' are 
time-integrals over the pulse envelope and its square. 

The population amplitudes a At) in Eq. ([2c|l can be ex- 
pressed as an integral of an(t) [26| and, after employing 
Eqs. (UJ), as an integral of a + (t) — a_(i). Using now the 
explicit expressions ([5]) makes the computation of a e (t) 
and of the spectrum cr(e) = |a e (oo)| 2 rather straightfor- 
ward 



a(e) = 



-r/4 j(t) 



|_ e ^ f+ Af( t )] +el [«-AF( f )]|^| 2 j (6) 



where we introduced the abbreviation S = IP + e — 2co — 
e — eq, which is the electron energy detuning from the 
center of the photoelectron spectrum e = 2w — IP. 

Interestingly, this expression for the spectrum can fur- 
ther be evaluated analytically. To this end we notice that 
the integrand in © contains the sum of two rapidly oscil- 
lating factors which is multiplied by a smoothly varying 
function of time. The main contributions to the inte- 
gral stem from the times at which two phases &±(t) = 
St±AF(t) are stationary i.e. $±(t s ) = 0. The two 
resulting stationary time conditions, 5 — ^fAg(t s ), have 
a transparent physical meaning. They define the time 
t s (e) at which an energy of a decoupled resonance, contin- 
uously shifted by the time-dependent coupling ±Ag(t), 
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FIG. 3: (Color online) Sequential two-photon ionization of 
H by a Gaussian-shaped pulse of 30 fs duration, carrier fre- 
quency of ui = 3/8 a.u. = 10.20 eV, and peak intensity of 
1.3 x 10 13 W/cm 2 . Shown are the photoelectron spectrum 
computed numerically (open circles; taken from Fig. [2]) and 
the spectrum obtained in the stationary phase approximation 
via the explicit expression Q (solid curve). The two individ- 
ual contributions to the spectrum, describing in (0 the sep- 
arate distributions of photoelectrons emitted at times when 
the pulse arrives and expires, are shown by broken curves. 



moves across the energy position 6 = e — Eq of the contin- 
uum state under inspection. These times and the elec- 
tron energy e are connected via the simple expression 
e = 2oj — IP =p Ag(t s ). During the pulse resonance |— ) 
covers the lower kinetic electron energy side of the spec- 
trum, e — £o G [ — A,0], and resonance |+) the higher 
energy side, e — £n € [0,+A]. For any pulse there are 
at least two stationary points for each value of e: one, 
ti(e), when the pulse is growing, and another, t2(s), when 
it decreases. For a Gaussian pulse there are exactly two 
times, ti(e) = -t 2 (e) = r v /ln[A/(e - e )]. 

By collecting in the integral (j6|) the two stationary 
phase contributions at t s — ±t±(e), we obtain the fol- 
lowing explicit approximate expression for the spectrum 



a(e) 



d£ 



-r/4j(t a 



t,=±ti(e) 



x |_ e i[# + (t»)Tf] +e i[*-(i s )±f]|| 2 . ( 7) 

The additional phase factors ? result from higher terms 
in the expansion of the phase <&±(t) around the station- 
ary points ±ti(e) computed for the Gaussian pulse [21 1. 
The photoelectron spectrum Eq. ([7]) is easily evaluated. 
The result is depicted in Fig. [3] by a solid curve. It is illu- 
minating to see that an explicit simple expression repro- 
duces nicely the numerically determined spectrum (open 
circles). The individual contributions of the two times 
t s = ±ti(e) to the spectrum in Eq. ([7]) are rather smooth 
and do not show any interference effects (broken curves). 

Eq. (O uncovers the physical origin of the strong mod- 
ulations in the electron spectrum. These are the results 
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FIG. 4: (Color online) The intense laser pulse of resonant 
carrier frequency induces a time-dependent coupling A g(t) 
between the ground and intermediate states. The energies of 
the resulting two decoupled resonances follow adiabatically 
the pulse envelope g(t) in two opposite directions (dashed 
curves). As a result, the photoelectron emitted by a second 
photon along the pulse envelope has at every moment t pre- 
dominantly a specific kinetic energy e. These are the times at 
which the energies of the decoupled resonances move across 
the energy position 5 = e — so of the continuum state. These 
passage times and the kinetic electron energy are simply con- 
nected via e = 2ui — IP =p Ag(t„). The pulse envelope first 
grows and then falls, and for a Gaussian pulse there are ex- 
actly two times at which the emitted electron wave has the 
same energy e. These two waves emitted with a time delay 
with respect to each other interfere, giving rise to the strongly 
modulated distribution of the photoelectrons shown on the en- 
ergy axis. Resonance |— ) is responsible for the low-energy part 
of the spectrum, and resonance +) for the high-energy part. 



of the coherent superposition of two photoelectron waves 
emitted with the same kinetic energy at two different 
times. A schematic visualization of the dynamic interfer- 
ence is given in Fig. 2) The dynamic interference spec- 
tacularly modifies the sequential two-photon ionization 
process and causes enormous qualitative changes in the 
spectrum, which can be verified by available laser pulse 
sources. The predicted effect is not constrained to se- 
quential two-photon ionization. We are convinced that 
dynamic interference is a very general and fundamental 
effect which is best manifested in the observable spectrum 
of the emitted particles by prominent multiple-peak pat- 
terns. Often, the dynamics of states coupled by intense 
laser pulses is governed by a Hamiltonian like that in 
Eq. ([3]) and this dynamics is in turn probed by emitted 
particles, either by employing an additional probe pulse, 
or by the same pulse. The emitted particles do not have 
to be photoelectrons. They can be, e.g., Auger electrons 
or photons. They all serve as a probe of the few-level 
system coupled by the pump pulse. 
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In the case of resonant Auger decay of an atom in a free 
electron laser field studied in Refs. [26|, |34j, |35[, a coher- 
ent high-energy laser pulse of resonant carrier frequency 
couples the ground state and a core-excited electronic 
state. The latter state decays by emitting an Auger elec- 
tron. The quantum motion of this two-level system is 
described by a 2 x 2 Hamiltonian similar to that consid- 
ered here (see Eq. (34) of Ref. 26}). The main differ- 
ence to the Hamiltonian ([3]) is the presence of a time- 
independent Auger rate Ta of the core-excited state in 
addition to the leakage by ionization to the continuum 
present in Eq. (J3j) - We thus have to add to — ^Tg 2 (t) de- 
scribing the leakage in Eq. ([3]) the term — a describing 
the Auger decay. The time-independent rate can be 
substantial 0, [13, [sH and is then not negligible com- 
pared to Ag(t) at least at the very beginning and very 
end of an intense pulse. However, whenever during the 
pulse the field-induced coupling Ag(t) between the two 
states becomes larger than the Auger decay width Ta, 
the above discussed scenario of decoupled resonances can 
be applied and dynamic interference takes place. Indeed, 
multiple-peak patterns in the Auger spectrum are found 
in Refs. [26l l34j . but hitherto not interpreted. In view 
of the present results, these patterns can be understood 
in terms of dynamic interference. In those cases where 
photons are emitted from the decoupled resonances, e.g., 
X rays, it is clear that the respective emission spectra of 
atoms exposed to coherent intense pulses will also exhibit 
dynamic interference effects. The only difference will be 
that a term — iTx will have to be added to Hamiltonian 
([3]) to account for the relaxation of the intermediate state 
via spontaneous emission. 
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